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m : ABSTRACT 



The magnetorotational instability (MRI) plays a key role for cosmic structure 
formation by triggering turbulence in the rotating flows of accretion disks that 
' would be otherwise hydrodynamically stable. In the limit of small magnetic 

Prandtl number, the helical and the azimuthal version of MRI are known to 
^ ■ be governed by a quite different scaling behaviour than the standard MRI with 

O ' a vertical applied magnetic field. Using the short-wavelength approximation 

^ ■ for an incompressible, resistive, and viscous rotating fluid we present a unified 

O ■ description of helical and azimuthal MRI, and we identify the universal character 

^ ' of the Liu limit 2(1 — \/2) ~ —0.8284 for the critical Rossby number. From 

■ this universal behaviour we are also lead to the prediction of higher azimuthal 

wavenumber for rather small ratios of azimuthal to axial applied fields. 

^ ■ Subject headings: instabilities, magnetohydrodynamics, waves 

.9. 

1. Introduction 

The magnetorotational instability (MRI) is widely accepted as the main source of tur- 
bulence and outward angular momentum transport that is needed for the matter in accretio n 



disks to spiral inwards onto the central proto-star or black hole, iBalbus fc Hawley I ( 119911 ). 
While the early work on MRI was mainly concerned with ideal magnetohydrodynamics, the 
last years have see n an increasing interes t in the influence of viscosity and electrical resis- 



tivity on the MRI, iPessah &: Chan I (|2008[ ). A particular role is thought to be played by the 
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magnetic Prandtl number (Pm) which measures the rat i o of v iscosity to resistivity. For ac- 
cretion disks around black holes (BH). iBalbus fc Henri I ( 120081 ) had discussed the transition 
from large values of Pm, in the vicinity of the BH, to small values in the outer part of the 
disk, with Pm reaching unity for approximately 100 Schwarzschild radii (depending on sev- 
eral parameters, among them the mass of the BH). By invoking a thermal runaway process 
at the unstable interface between regions with Pm > 1 and Pm < 1, the authors associated 



this b oundary with the existence of high and low X-ray states, see iRemillard fc McClintock 
J2006h . 



In recent yea r s, a v i vid discussion (ILesur &: Longaretti I (120071 ) ; iFromang et al. I ( 120071 ) ; 



Kapyla fc Korpi I (l201ll ): lOishi &: Mac Low I (l201ll )) was devoted to the possible decline of 
the angular momentum transport rate with decreasing Pm, and to the intricate roles that 
are played here by the magnetic Reynolds number, the detailed boundary conditions, and 
the stratification of the disk. 

Besides this relevance to the outer part of accretion d isks, to protoplaneta r y dis ks 
(iTurner fc Sano I (120081 )). and possibly even to planetary cores ( jPetitdemange et al. I ( l2008l )). 
the limit of low Pm has acquired some additional i nterest in connection with the recen t liqui d 



metal experimeii t s devo ted to the study of MRI (ISisan et al. I (120041 ): IStefani et al. I ( 120061 ): 



Nornberg et al. I ( 120101 )). While for the standard version of MRI (SMRI), characterized by 
only a vertical field bei ng applied, the low Pm limit is rather smooth and unspectacular 
( Pessah &: Chan ( 20081 )). the addition of an azira u thal fi eld leads to dramatic effects as re- 
vealed for the first time by iHoUerbach &: Riidiger I (120051 ). The arising helical MRI (HMRI), 
as we call it now, was shown to work also in the inductionless limit since it is governed by 
the Reynolds and Hartmann number, quite in contrast to SMRI which is governed by the 
magnetic Reynol ds number and the Lundquist number. However, as it was early shown by 



Liu et al. I ( l2006l ). the functioning of HMRI is limited to comparably steep rotation pro- 
files with Rossby numbers Ro < ^ouu = 2(1 — \/2) ~ —0.8284 (which we henceforth will 
call the "Liu limit"), and does therefore not extend to the astrophysically important Ke- 
pler profiles characterized by Ro = —0.75. This essential limitation of the HMRI, together 
with a va riety of further parameter d epend encies, was confirmed in the PROMISE exper- 



iment by IStefani et al. I (120061 120071 . l2009l ). The intricate, though continuous, transition 
between SMRI and HMRI which involves a spectral exceptional point at which the inertial 
wave branc h coalesces with the bran c h of t he slow magnetocoriolis wave, was clarified only 



recently by iKirillov k Stefani I (120101 120111 ). 



Another surprise in the limit of low Pm was, for the case of a purely or strongly dom- 
inant azimuthal magnetic field, the numerical prediction of a no n-axisymmetric version of 
MRI, working apparently in a similar parameter region as HMRI ( IHoUerbach et al. I ( 12010! )). 
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Although the occurrence of MRI ca se under the influence of a pure l y azimuthal magnetic field 
had b e en studied much earlier, see iHawlev. Gammie and Balbus I ( 119951 ) ; lOgilvie &: Pringle 
( 119961 ) ; iTerquem &: Papaloizou I ( 119971 ) , the crucial effect arose again for the particular com- 
bination of low Pm and slightly steeper than Keplerian shear profiles. It has to be noticed 
that this azimuthal MRI (AMRI), as we call it now, works for azimuthal magnetic fie lds that 
are cu r rent-free in the consider ed fiuid, quite in contrast to the Tayler instability (ITayler 
( 119731 ): ISeilmayer et al. I (120121 )) that is a pinch- type instability in a current-carrying con- 
ducting medium. 

The aim of this paper is to better understand why the scaling behaviour of HMRI 
and AMRI, and their restriction to rather steep rotation profiles, is so similar. In order 
to clarify this point we restrict our work completely to the short wav elength approximation 
(IFriedlander &: Lipton-Lifschitz I ( 120031 ): iHattori fc Fukumoto I (120031 ) ). keeping in mind that 
some of our conclusions will need further confirmation in more realistic simulations. 



2. Short wavelength analysis of viscous, resistive MRI for arbitrary azimuthal 

wavenumbers 

We start from the equations of incompressible, viscous and resistive magnetohydrody- 
namics, comprising the Navier-Stokes equation for the velocity field u and the induction 
equation for the magnetic field B, 

BVB VP ^2 

+ z/V^u (1) 

Pop P 

B- Vu-u- VB + r^V^B, (2) 

where P = p + ^ is the total pressure, p = const the density, u = const the kinematic 
viscosity, r] = (/ioci)"^ the magnetic diffusivity, a the conductivity of the fiuid, and po 
the magnetic permeability of free space. Additionally, the mass continuity equation for 
incompressible flows and the divergence-free condition for the magnetic induction are used: 

V-u = 0, V-B = 0. (3) 

In the following we consider a rotational flow in the gap between the radii Ri and 
i?2 > -Ri, with an imposed magnetic fleld sustained by currents external to the fluid (hence 
we disregard any version of the Tayler instability and its combination with MRI). Introducing 
the cylindrical coordinates {R, 0, z) we consider the stability of a magnetized Taylor- Couette 
(TC) flow, i.e. a steady-state background flow with the angular velocity proflle fl{R) in a 



- + U.VU 

dB 

'dt 
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(generally helical) background magnetic field, 

uo{R) = Rn{R) e^, p = Po{R), Bo(i?) = 5°(i?)e^ + 5°e, 
with the azimuthal field component 



B'AR) 



2tiR' 



(4) 



(5) 



supposed to be produced by an axial current I. 

The angular velocity profile of the background TC fiow is known to have the form 



Q{R) = a+—, 



where a and b are defined as 



a = — TT^isi, 



1 - p2 

1 — 17^ 



Ri^i 



(6) 



(7) 



with the definitions 



Ri 



^2 



Introducing, as a measure of the steepness of the rotation profile, the Rossby number (Ro) 

R dn 



we find 



2Q OR 

a = fi(l + Ro), h=-VLR^Ro. 



(9) 
(10) 



To study fiow and magnetic field perturbations on the background of the magnetized TC 
fiow we linearize the Navier-Stokes and induction equation in the vicinity of the stationary 
solution by assuming u = uq + u', p = + p', and B = Bq + B' and leaving only terms of 
first order with respect to the primed quantities. 

Then, by using a short-wavelength approximation (the details of the d erivation will be 



publis h ed elsewhere) in the frame of the geometrical optics approach (see e.g. iLandman fc Saffman 



1987 ): 



fl2003h : 



Dobrokhotov fc Shafarevich l(11992r):lFriedlander fc Lipton-Lifschitz I (120031 ) : lHattori fc Fukumoto 



Lebowitz fc Zweibel I (120041 ): iMizerski fc Bajer I ( 12 0091 )) we end up with a system of 4 



coupled equations for the perturbations of arbitrary azimuthal dep endency which generalize 
the corresponding equations derived in iKirillov fc Stefani I (120101 ). From those 4 coupled 
equations, we can deduce the dispersion relation 



p(7) := det(i/ -7^) = 



(11) 
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generated by the matrix 
/ 



H 



\ 



—imVt — LOi, 
-2Q{1 + Ro) 



2a^n 
—irriQ — 








—imVt 



UJ„ 







121 



where we have used the following definitions for the viscous, resistive, and the two Alfven 
frequencies corresponding to the vertical and the azimuthal magnetic field: 



Z/|k|2 
|2 



^1 



?7|k| 
klB 



02 



(13) 
(14) 

(15) 
(16) 



Note that |kp = kj^ + k^, and a = kz/\^\, where k^, m, and kz are the radial, azimuthal, and 
axial wavenumbers of the perturbation. In the absence of the magnetic field , the dispersion 
relatio n determined by the matrix H reduces to that derived already by iKrueger et al. 
(1l966l ) for the non-axisymmetric perturbations of the hydr odynamical Couet t e- Tay lor fiow. 
Choosing, additionally, m = 0, we reproduce the result of lEckhardt &: Yao I (119951 ) . In the 
presence of the magnetic fields and m = 0, we arrive at the dispersion relation derived by 
Kirillov fc Stefani I fcoiol ). 



The dispersion relation flTTl) generated by the matrix f[T^ can be rewritten completely in 
terms of dimensionless numbers, i.e. Rossby number (Ro), magnetic Prandtl number (Pm), 
ratio of the two Alfven frequencies {(3), Hartmann number (Ha), Reynolds number (Re) and 
a re-scaled azimuthal wavenumber n: 



(17) 
(18) 
(19) 

(20) 
(21) 



Pm = 






V 


/3 = 




a — - 






n 


Re = 


a — 







Ha 



5 



a 



ky/lIopUT] 



n = — 



m 
a 
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After re-scaling the spectral parameter as 7 = X^/uJi/j^ we end up with the complex 
polynomial dispersion relation 

p{\) = aoA^ + (ai + ibi)\^ + (aa + ib2)X'^ + (03 + ih)>^ + ai + ih = (22) 

with the coefficients: 

ao = 1 

V Vp^J 

bi = 4nRevPm, 

02 = 2(/32Ha2-3Re2Pm)n2 + 4/3Ha2n + 2(l + (l + 2/32)Ha2)+4Re2(l + Ro)Pm + ^ 

62 = 6nRe(l + Pm) 

ag = ai(/32Ha2-3Re^Pm)n2 + 2ai/3Ha^n + ai(l + (l + 2/32)Ha2) + 8Re2(l + Ro)yP^ 

63 = 4n3yP^Re(/32Ha2 - Re^Pm) 

+ 2nRe(4Pm2Re^(l + Ro) + (1 + Pm)^ + 2Pm(l + Ha^))/^?^ 

m 

04 = ((/3W - Re^Pm)^^ + 2Ha2/3n + Ha^ + 2PmRe^)2 

+ 2(2Re^PmRo + l){(Ra^l3^ - Re^Pm)n2 + 2Ha2/3n + Ha^) - (1 + FmfRe'^n^ 

+ 4Re2(l + Ro) - (Ha^ + 2PmRe2)2 + Ha^ + 1 + AjS^Ra^ 

64 = 2Re(l + Pm)(/3W-Re2pm)n^ + 4ReHa2/3(l + Pm)n2 

+ 2Re(2(l + Ro)(2Re2pm - l3^}ia^{l - Pm)) + (1 + Ha2)(l + Pm))n 

- 4/3Ha^Re(2+ (l-Pm)Ro). (23) 



Note again that this complex algebraic equation of 4th order is valid for perturbations 
of arbitrary azimuthal wavenumber in magnetized incompressible, viscous, resistive rotating 
fluids exposed to current free axial and azimu thal magnetic fields . Wh en n = it reduces 
to the dispersion relation of HMRI derived by iKirillov fc Stefani I (120 lOf ). 



Inductionless limit 



Proc e eding quite similar as in iKirillov fc Stefani I ( I2OIOI ) , we apply the Bilharz criterion 
( iBilharz I (jl944( )) to the complex polynomial ( 122|) and derive the maximum Rossby number, 
at which flows are prone to MRI, as a function of the remaining dimensionless numbers. In 
the following, we concentrate on the inductionless limit, i.e. we take the limit Pm — )■ 0. 



Fig. 1. — Inductionless limit Pm — )■ 0: (a) At the given n = y/2 - 1/20 and (3 = 20 the 
threshold of MRI in the (Re~^, Ha~^, Ro) space with the maximum at the singular (red) 
point (0,0,Rocr = 2 — 2v^). The red line leading to the maximum projects into a curve 
approximated in the vicinity of the origin by the scaling law ( 125|1 . (b) The threshold of MRI 
in the Ro) space with Ha = 6 and Re determined by the scaling law ( 125|) reaches its 
upper bound at /9 — )• oo (AMRI), which is still below the Liu limit Rouu = 2-2^2. 



After verifying the facts that the inductionless threshold value of Ro increases monotoni- 
cally with Re so that we can take the limit Re — )■ oo, and that in this limit the threshold value 
of Ro increases monotonically with Ha, Fig. [D^a), so that we can take the limit Ha — )■ oo, 
we obtain the following explicit expression for this maximized (with respect to Ha and Re) 
critical Rossby number at the threshold of MRI in the inductionless limit: 



(o . _ 4/3^ + {pn + ir- {2/3^ + {/3n + l)') ^4/3^ + (/3rz + 1)^ 
t^Ocr[P,n) - 2/32(/3n + l)2 ' 

The maximum value of the Rossby number, Rocr, is shown as a red dot in Fig. [D^a). When 
n = 0, Eg. f[2^ r educe s to the threshold of HMRI in the inductionless limit found in 



Kirillov fc Stefani I (120111 ). Given n and /3, the critical Rossby number calculated with fi- 
nite values of Re and Ha, is always below the majorating value, Rocr-(/3, '^), determined by 
Eq. (^^, see Fig. [T](b). With the increase of Ha and Re constrained by the scaling law 

Re = 2(1 + y2)/3=^Ha^ (25) 

the threshold shown in Fig.[^b) tends to the majorating surface Rocr(/3, n) shown in Fig.[2](a). 



The remarkably simple dependence fl24l) of Rocr, only on the ratio (3 of azimuthal to 
axial field and on the rescaled azimuthal wavenumber n, relies on the appropriate choice of 
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the dimensionless parameters, in particular on "hiding" the wavenumber ratio a in them. 
Assume we have fixed the sign of /3, the functional dependence of the maximum critical 
Rossby number f l2^ on /3 and n has a two-saddle line structure as shown in Fig. M^a). 




b) 

4 



-2 



















■ T 1 1 1 1 1 1 1 1 1 — 


i ( P 



-2-1012 



Fig. 2. — (a) The critical Rossby number maximized with respect to Ha and Re, given by 
Eq. flMl) . in dependence on /3 and n. (b) The lines (green) P>raax,\ and (red) (3max,2 at which 
the function KocriP^n) attains its maximal value Ro^j^ = 2 — 2\/2. 



Now, the most important result of this paper is that both saddle lines have the same 
height everywhere, namely the Liu limit HoLiu = 2(1 — a/2), independent on the particular 
combination of (3 and n. From Eq. (l24l) it can be proved that these saddle lines, with 
RoLiu = 2 — 2\/2, are governed by the two equations 



Pmax,l ('^) 
I3max,2{n) 



—= , ROcr(/3maa',l) = 2 — 2a/2 

V 2 — n 



-1 



72 + 



n 



ROcr(/3. 



maxfl I 



2 - 2V2. 



(26) 



In Fig.El^b) the curves [imax^iip^ and l3max,2{^) are shown in green and red colors, respectively. 
Note that according to the first of Eqs. fl26|) . n = (HMRI) corresponds to /3 = 1/a/2 which 
being substituted into E g. (I25l). yields the fo l lowing scaling law for the optimum combination 
of Re and Ha in HMRI kirillov fc Stefani I kom ) 



Re 



(27) 



Therefore, even in the case of non-axisymmetric perturbations, the maximum possible 
value of the Rossby number prone to the magnetorotational instability caused by the helical 
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magnetic field in the inductionless limit is still HoLiu = 2 — 2\/2, exactly as in the case of 
HMRI which is an instability with respect to the axisymmetric perturbations (n = 0). The 
relations (l26l) between /3 and n that correspond to the Liu limit give a sort of the resonance 
conditions between the components of the wavevector of the three-dimensional perturbation 
and the components of the helical magnetic field. 

1 

0.8 
0.6 
0.4 
0.2 


0.5 1 1.5 2 2.5 3 
P 




Fig. 3. — Range of a for the different m modes, in dependence on /3. 

On the basis of the equations ( 126|) connecting /3 and the rescaled azimuthal wavenumber 
n, we can ask now for the structure of the solution in terms of the original, "physical" 
azimuthal wavenumber m. From the definition a = A;^/|k| we see immediately that a can 
take on only values between -1 and +1. The solution structure can thus be visualized as in 
Fig. [31 For example, from the first of the Eqs. we see that for large values of (3 (AMRI), 
the only possible integer solution is the m = 1 mode, whose corresponding wavenumber ratio 
converges than to a = see also Fig.[2t^b). There is a lower limit of /3 = 1+^2 ~ 2.41 for 

this m = 1 mode. Lowering (3 further, we find n ext the m = mod e (HM RI) to dominate 



at the Liu limit, restricted only to /3 = 1/v^, Kirillov fc Stefani I ( 20101 ). Interestingly, 
decreasing /3 > even further to zero we find a sequence of higher azimuthal modes with the 
sign of m that is opposite to the sign of (3, which indicates a kind of resonance phenomenon. 

Since a enters also the definition of f3 it might be instructive to illustrate the mode 
structure also in dependence on (3/ a = coa^/'^a- From Eqs. (12^ we derive 

a = ±^L + ^), (28) 

where the positive sign corresponds to the first of Eqs. f l26|) and the negative sign to the 
second one. This means that with a given azimuthal wavenumber m, two axial wavenumbers. 
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kz, are associated following from Eq. fl28|l that differ by sign only. Such combinations of 
wavenumbers are the most destabilizing in the sense that the magnetized Taylor- Couette 
flow is unstable at the highest possible Rossby number. 

In Fig. m we plot the positive branch of Eq. (128|) because the negative one is simply its 
reflection about the horizontal coordinate axis. We see now the HMRI mode (m = 0) to 
start at uja^/(^a = 1/v^ and to remain for arbitrary large values of (3/ a, although with an 
ever decreasing wavenumber ratio, which would correspond to ever increasing wavelengths 
in z direction. Again it is only the AMRI modes (m = ±1) that, for large P/a, maintain a 
physically sensible wavenumbers a = ±1/a/2. The higher modes with m < —2 are obtained 
for smaller values of 



-1 





1 \ ^ 


J 


\-2 

-3\ ' 



2 4 



Fig. 4. — Range of a for the different m modes, in dependence on (red line) m = 1, 

(black line) m = 0, (green lines) m = —1, —2, —3, —4, —5. 

Hence, when the azimuthal magnetic field is directed along the basic flow that rotates 
counter-clockwise with respect to 2;-axis, among the modes that are MRI-unstable at the Liu 
limit there are only two (AMRI) that co- rotate with the flow (m = 1) and simultaneously 
propagate either along the positive or negative z-direction. These modes are dominant when 
CO /to A > 1 + V^, see the red curve in Fig. HI At moderate ratios to a^/ to a, the Liu limit is 
at the axisymmetric HMRI mode. When coa^/^oa — ?■ 0, infinitely many modes with m < —1 
that propagate either in the negative or positive z-direction and counter-rotate with respect 
to the basic flow can cause instability at the Liu limit. Note, however, that at finite Re and 
Ha the highest modes will be inhibited, see Fig. [I](b). 
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Conclusion 



Using a sliort-wavelengtli approacli, we liave presented a unifying picture of the induc- 
tionless forms of MRI. We have identified a continuous function of the maximum critical 
Rossby number that incorporates both types of instabihty. We were lead to the conclusion 
that in the limit of small ratios of azimuthal to axial field there should be inductionless MRI 
versions with higher m modes, counter-rotating with respect to the basic flow, although this 
needs further confirmation at least by a ID linear stability analysis. Most interestingly, the 
Liu limit has turned out as being of quite universal significance, since the range of its validity 
has been extended from the realm of axisymmetric HMRl to that of non-axisymmetric MRI 
versions. Actually, soon after its derivation in the WKB approxi mation, the relevance of the 
Liu limit had been questioned by lRiidiger and HoUerbach I ( 120071 ) who had found an apparent 
extension of this limit in global simulations when at least one of the radial boundary condi- 
tions was assumed to be electrically conducting. Later, though, utilizing another definition 
of the notion "quasi-Keplerian" for Taylor-Couette flows, the Liu limit was rehabilitated by 
Priede I (l201l[ ). As a side remark, the determination of such limits is even more complicated 
by the necessity to distinguish, for travelling waves as in HMRl, betw een convective and 



absolu te (or global) instabilities, which has been thoroughly disc ussed bvlPriede fc Gerbeth 
( 2009 ) and which was shown to be experimentally important by Stefani et al. ( 20091 ). 



From the strictly astrophysical point of view, our support for the Liu limit may appear 
disappointing, since it would exclude any relevance of the inductionless versions of MRI 
to accretion disks with Keplerian rotation. A subtly question in this respect is, however, 
connected with the saturation mechanism of the MRI that could , possibly, lead to modified 
shear profiles. With main focus on low Pm flows, lUmurhan I (120101 ) had asked for the 
possibility that the saturation of MRI could lead to modified flow structures within parts 
of steeper shear, sandwiched with parts of shallower shear. By virtue of a possible sudden 
onset within such segments of steepening shear, the inductionless MRI versions could thus 
play a certain role in real astrophysical settings. 
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